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Abstract 


Thermal structure calculations have been carried out for the 
atmosphere of Jupiter above the ^evel of the dense clouds, for boundary 
conditions suggested by recent observations* The resulting models are 
characterized by an extensive region of dynamical control above the 
cloud level, and a thermal inversion in the mesosphere, produced by 
absorption of solar infrared energy in the 3020 cm~^ band of methane. 

The infrared and microwave spectra corresponding to the computed thermal 
models are found to be in generally good agreement with observed infra- 
red and microwave brightness temperatures of Jupiter, The effective 
temperatures of all of the thermal models are higher than the solar 
eouilibrium value, and, thus, an internal heat source on Jupiter is 
suggested , 


Introduction, 


During the past year, important new observations of Jupiter have 
provided the motivation for a reinvestigation of the thermal structure 
of the Jovian atmosphere. First, several authors have recently reeval- 
uated the abundances of hydrogen, methane and ammonia in the atmosphere 
of Jupiter, by estimating the total amounts of these constituents above 
various temperature and pressure levels (Owen and Mason, 196$; Belton, 
1969). It now appears that the atmospheric composition on Jupiter is 
similar to that of the sun, and that the hydrogen— helium ratio by 
volume may be as large as 5.0. Second, the brightness temperature of 
Jupiter has recently been measured at a number of infrared wavelengths 
between 2 ^ and 20^ by Gillett, et al. (1969). Important suggestions 
of those observations are (a) that the brightness temperature of Jupiter 
in the atmospheric window near 5 ft may be as high as 230K, and may be 
characteristic of the temperature at the top of a thick cloud, and (b) 
that a temperature inversion may be present in the middle atmosphere. 
Finally, microwave measurements of Jupiter have been carried out at 
several wavelengths in the vicinity of the strong absorption by the 
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numerous inversion lines of ammonia in the 1 cm to 2 cm region (law and 
Staelin, 1966; Welch, 1969). These measurements in the center and wings 
of the NH-j inversion band provide further information on the atmospheric 
temperature distribution. 

These new sets of data constitute separate constraints which must 
be satisfied by any acceptable model of the vertical structure of the 
Jovian atmosphere. A theoretical model should be consistent with the 
latest estimates of atmospheric composition, while leading to a spectral 
distribution of radiative loss to space in agreement with the observed 
infrared and microwave brightness temperatures. 


We wish to renort here on a calculation of temperature and density 
profiles for the atmosphere of Jupiter above the level of the dense 
clouds. Thermal structure calculations have been made for various com- 
binations of temperature, total pressure and composition at this level. 
For each thermal model, the energy flux radiated by the planet at 
various infrared and microwave wavelengths was also computed and com- 
pared with observations. 

Vertical Structure . 

In the present study, a steady-state distribution of temperature 
and density for the atmosphere of Jupiter above the cloud level is com- 
puted as follows. Basic atmospheric properties at the cloud level, 
including the cloudtop temperature, the mixing ratios of the major 
constituents (F 2 , He, CH^ t NH^) and the total pressure are adopted, and 
a zero-order vertical temperature profile above the clouds is specified, 
A consistent distribution of pressure and density is then determined, 
assuming the atmosphere to be in hydrostatic equilibrium. The concen- 
trations of H 2 , He and CH^ are computed with the assumption of a mixed 
atmosphere, and the mixing ratio of NH^ is assumed to be constant at 
levels where NH^ is unsaturated. At levels where NH^ saturation occurs, 
the concentration of NH^ is computed from a saturation formula derived 
from Lasker (1963) 

”nh (t) * 

3 exp [ 25.88 x . 1.0)1 

where T is the temperature and n (145.0) is the saturation number 


(145.0) x ( ) 


density of NH^ at 145K, taken here to be 2. #4 x lO 1 ^ cm“^. 

For the zero-order structure obtained in this manner, the mean 

infrared transmission of the atmosphere in wavelength intervals of 

interest is evaluated for required paths, considering contributions 

to the infrared opacity by pressure- induced absorption of H 2 and the 

vibration -rotation bands of CH and NH 0 . Using these transmissions, 
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the energies of planetary - 100^) and solar (2^ - 4^) origin 
absorbed per unit volume at various levels in the atmosphere are com- 
puted, and the total energy absorbed is equated to the ^ - 100^ 
emission. The energy balance relationship is thus 

e,z) + S(v,e,z)j ^n i (z)k i (v> > z) d»d« 

* JfJjC^.zlT 1 . n^( z)k^(v, z) dv>dco 


where l(y,0,z) is the intensity of planetary radiation, S(y,0,z) is the 
solar intensity, J(v,z) is the atmospheric source function, n^(z) and 
k^(v,z) are the number density and absorption coefficient of the i th 
constituent, and V, ©, z and co have their usual meanings. 

In the above expression, I(y,©,z) and S(V,0,z) are evaluated by 
standard radiative transfer techniques. It is assumed that the cloudtqj 
radiates like a blackbody in the 7// - lOOytf range, and that the downward 
directed intensity of planetary radiation at the upper boundary is 
negligible. S(y,0,z) is computed for daily mean conditions at latitude 
45®, equinox. 

By relating the 7 /{ to lOO^f emission (represented by the right hand 
side of the above equation) to the temperature, assuming that local 
thermodynamic equilibrium conditions prevail at all levels under con- 



sideration, a first order profile of temperature is obtained. The 
atmospheric structure is then recomputed for the first-order temperature 
distribution. Iteration of this procedure is performed until tempera- 
ture convergence to within an arbitrarily small limit is achieved. 

In the course of these calculations , whenever a superadiabatic 
temperature gradient is obtained , the calculated lapse rate is replaced 
by the dry-adiabatic lapse rate. In this manner, allowances are made 
for adjustments of unstable radiative profiles by dynamical processes in 
the optically thick region of the Jovian atmosphere adjacent to the 
clouds. 

The release of latent heat in the NH^ sublimation process would nob 
significantly alter the basic thermal structure of the atmosphere, and 
this effect was not taken into account in the energy balance relation- 
ship. 

In the present solution of the equation of radiative transfer, the 
infrared regions between 2 and 4 f\ and between and 100 J 4 were 
divided respectively into sixteen and fifty wavelength intervals of 
Varying width. Higher resolution was adopted in regions where rapid 
variations of absorption coefficient with wavelength occur, tc obtain an 
acceptable simulation of the contours of absorption both in the CH^ and 
NH^ vibration-rotation bands and in the continuum. 

Kean absorption coefficients for CH^ in the 3020 cnf^ and 1306 cm"^ 
bands were estimated from the measurements of Burch and Williams (1962), 
while absorption coefficients in the 3300 cm“^, 1628 cm“^ and 950 cm“^ 
bands of NH^ were estimated from the measurements of France and Williams 
(1966). In extrapolating the CH^ and NH^ absorption coefficients to 


obtain their values under atmospheric conditions, the strong line 
approximation (Plass, i960) was employed. 

Coefficients for Hg absorption, induced under high pressures by 
H 2 - H 2 and Hg - He collisions, were taken from Trafton (1967). 

Traftcm’s values were extrapolated to atmospheric conditions taking 
into aoeount the strong pressure dependence (pressure 2 ) of absorption 
of this type. 

Values of the solar radiation in the 2j\ - 4yl< region incident on 
the top of the earth’s atmosphere were obtained from the Handbook of 
Geophysics (Gast, I960). In correcting these values for the orbital 
distance of Jupiter, the mean orbital radius of 5.203 A.U. was adopted. 

Five thermal models of the Jovian atmosphere were constructed by 
the method outlined above, for cloudtop temperatures ranging from 210K 
to 230K, and for total pressures at the cloud level between 2.4 atm and 
7.2 atm. The above range of temperature at the cloudtop is consistent 
with the findings of Gillett, et al. (1969), and also with the results 
of Owen (1965) and Danielson (1966), who deduced temperatures of 
200 t 25K and 210 - 15K, respectively, in the vicinity of the cloudtop. 
The proportions by volume of H2, He, CH^ and NH^ were assumed to be 
5.0 : 1.0 : 0.005 t 0.001. Basic characteristics of these model 
atmospheres are listed in Table I, and the temperature profiles obtained 
in these cases are shown in Figure 1. 

In all five of the resulting thermal structures, a rather extensive 
region of dynamical control develops above the cloud level, with an 
adiabatic decrease in temperature from the cloudtop to a "tropopause" 
in the 30 km - 40 km altitude range. Above the "tropopause" , tempera- 
tures continue to fall sub-adiabatieally to a low value near the 65 km 



level, but in all cases except Model 1, an increase in temperature with 

altitude is found at levels above 65 km, as a result of near infrared 

solar heating in the 3020 cm CH band. As one might anticipate, the 
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strength of this inversion is found to depend directly, on the cloudtop 
pressure (opacity) for a given cloudtop temperature* and inversely on 
the cloudtop temperature for a given cloudtop pressure (opacity). No 
inversion appears in Model 1, but, rather, a general warming of the 
atmosphere at levels above the "tropopause" is found* as a consequence 
of the relatively low pressure, small fai 1 'rared opacities, and warm 
radiative equilibrium temperatures which characterize this model. 

Gillett* et al* (1969) interpreted their infrared observations of 
Jupiter as indicating the presence of an inversion in the Jovian atmos- 
phere* arid Were the first to suggest that solar heating in the 3020 cnf* 
CH^ bahd could be responsible for its formation. Thus* the present cal- 
culations confirm their hypothesis. Solar Ultraviolet photoionization 
and photodisSOciation energy is deposited well abcve the 100 km level in 
ail of the Models considered here* and does not aff.ect the thermal 
structure of the atmosphere in the altitude range of interest* 

The vertical distribution of Hg, tie, and NH^ corresponding to a 
typical model (Model 3) are Shown in Figure 2. DUe to their extreme 
similarity to the pattern of ?igure 2,' the density distributions for the 
other thermal models have been omitted. The sharp decrease in the NH3 
concentration betweeh 30 km and 45 km altitude is produced by NH-j satu- 
ration near tha 145K level. The maximum amount 6f Nti^ Which could be 
deposited in this regiCn is O.65 gm cm" 2 * the difference between the 
amount which would be present above the level of saturation if the sub* 
limation process did not occur, and the amount present on the basis of a 


saturated NH^ distribution. Even if all of this crystalized NH^ were 
distributed uniformly over a 15 km altitude range, the correspond inr 
density of solid material would be only 4.0xl0~^ gm cm - ^ , of approxi- 
mately the same magnitude as the density of cirrus clouds in the earth’s 
atmosphere. 

Similar estimates of the density of a thin cloud of solid NH^ on 
Jupiter have been obtained by Lewis (1969). It appears unlikely that 
the presence of such clouds would modulate the transfer of far infrared 
radiation to any large extent. Indeed, the high brightness temperature 
observed near 5/< by Gillett, et al. (1969) implies that the transmit- 
tance of these clouds is quite high at this wavelength. It is not un- 
reasonable, therefore, to assume that they are equally transparent at 
longer wavelengths. 

In Table II, the total amounts of Hg, He, CH^ and NHj present above 
various temperature levels in the atmosphere of Model 3 are listed. 

These abundances agree in general with the estimates of Savage and 
Danielson (1968), Owen and Mason (1968) and Belton (1969). 

The effective temperature of the radiative losses to space corres- 
ponding to thermal Models 1 through 5 are 145K, 132K, 127K, 121K and 
112K, respectively. A value of 140K for the effective temperature of 
Jupiter was measured by Low (1969). Although the effective temperatures 
of the present models are generally somewhat lower than that observed by 
Low, they are, nevertheless, higher than the range of 88K to 105K, the 
solar equilibrium effective temperature corresponding to the range in 
the measured albedo of Jupiter of 0,73 (Harris, 1961) to 0.45 (Taylor, 
1965). Our models, therefore, imply the existence of an internal heat 
source on Jupiter. 


Atmospheric dynamics were taken into account in the present study 
only insofar as they might adjust the radiative structure of the atmos- 
phere in unstable regions, via convection. No rigorous treatment of 
Jovian dynamics was included in the one-dimensional, steady-state prob- 
lem considered here. However, the fact that an extensive convective 
zone develops in all of the present models, allows certain conclusions 
to be drawn about the importance of radiative transfer in determining 
the structure of the dense regions of the Jovian atmosphere adjacent to 
the clouds. From these calculations, it is clear that radiative trans- 
fer, acting alone, would produce a deep region characterized by a supers 
adiabatic temperature gradient. It is thus indicated that radiative 
coupling of atmospheric layers must be weak for a considerable distance 
above the cloud level. It is likely, therefore, that dynamical proces- 
ses of one type or another would strongly influence the thermal struc- 
ture in this region. 

From their analysis of the relative magnitudes of the time con- 
stants for radiative and dynamical processes, Gierasch and Goody ( 1969 ) 
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also concluded that the structure of the Jovian atmosphere is determined 
mainly by dynamics over an extensive altitude range. Their findings, 
which are based upon a more thorough analysis and direct approach to the 
problem, are supported by the results of the present study. 

Regardless of the nature of the dynamical processes which determine 
the structure of the atmosphere above the cloud level, the end result 
must be the production of a zone of adiabatic temperature gradient, sim- 
ilar to the "troposphere" which develops in the present models. The 
thickness of the "troposphere" would depend upon the strength of the 
circulation responsible for its formation, but could be no less than the 
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thickness computed here, on the basis of a convective adjustment of the 
radiative temperature profile. Cn the other hand, were this zone of 
dynamical influence to extend to much higher altitudes, the atmospheric 
inversion would be suppressed by dynamical redistribution of the thermal 
energy. In that case, no atmospheric inversion would be in evidence in 
the infrared spectrum of Jupiter. Therefore, it appears unlikely that 
dynamical processes would greatly modify the structure of the Jovian 
atmosphere from the basic pattern of Figure 1. 

Infrared Brightness Temperature . 

For each of the five thermal models of the Jovian atmosphere, the 
radiative flux emanating from the planet at various infrared wavelengths 
was calculated from radiative transfer theory, and a brightness temper- 
ature at each wavelength was derived. In the determination of bright- 
ness temperature, the transmission of the atmosphere at wavelengths of 
interest was obtained in the manner previously described. 

In Figure 3, the resulting infrared brightness temperatures for the 

* 

five thermal models are shown. The basic patterns which emerge are sim- 
ilar in all five cases, and in those spectral regions where near over- 
lapping of several brightness temperature curves occurs, some of the 
curves have been omitted from the figure to avoid confusion. 

The low value of brightness temperature found between 3^ and 4 ^ 
is a result of strong absorption in the 3020 cm"* CH^ band, while the 
smaller feature at 3/f is due to the 3300 cm"* band of NH^. The 162# 
cm"* NHj band produces the lew brightness temperatures near while 

the 1306 cm“* CH^ band provides the opacity in the 7/< - 6,5 ^ region. 
Atmospheric absorption between £.5^f and 12 is mainly due to the broad 


950 cm ^ band of NH^, but includes a component due to Hg which increases 
with wavelength in this range. At wavelengths greater than 12 A , the 
NHj opacity decreases rapidly, and Hg absorption is primarily respon- 
sible for the computed brightness temperatures beyond 12 yq • The H2 
opacity reaches a maximum near 17 /* , and decreases again at longer 
wavelengths . 

The effect of the atmospheric inversion on the calculated bright- 
ness temperatures can be discerned in three spectral regions where 
strong absorption is present, near the centers of the 3020 cm”'*' CH^ 
band, the 162$ cm ^ NH^ band, and the 1306 cm ^ CH^ band. In the strong 
regions of these bands, the highest brightness temperatures are found in 
the most opaque models. In Model 1, a "low opacity" model where no in- 
version is present, local minima in brightness temperature are computed 
at 3.3/f and 7*7/1 near the centers of the CH^ bands. However, in Model 
3, for example, a model of medium opacity where an inversion develops, 
local maxima are found at these wavelengths. 

Also plotted in Figure 3 are points corresponding to the measure- 
ments of Gillett, et al. (1969) and Low (1969). Earlier observations of 
brightness temperature in the - 14 /| range have resulted in values 
between 120K and 140K (Menzel, et al., 1926; Sinton and Strong, I960; 
Murray and Wildey, 1963; Sinton, 1964). As. can be seen, the calculated 
values are in general accord with the observations in the 7/< - 14 yt< 
range, with a tendency toward the observed elevation in brightness tem- 
perature between 7yq and 9/< due to the atmospheric inversion in Models 
2 through 5* Good agreement with the measured brightness temperatures 
in the atmospheric "window" between UJi and 5/f is also obtained. Best 
overall agreement with the observations is found in the case of Model 3* 


As mentioned earlier, in all of the models presented here, the same 
atmospheric composition was adopted (Hg : He : CH^ : NHj - 5.0 : 1.0 : 
0,005 : 0 . 001 ), Agreement with the observations of Gillett, et al. 
(1969) could undoubtedly be improved by adjusting the relative concen- 
trations of these constituents. For examole, a change in the concentra- 
tion of NH^ at higher levels in the models, perhaps due to the use of an 
alternate NH^ saturation relationship, would modify the calculated 
pattern of brightness temperature in the #. 5 /* - I 2 y^ wavelength region, 
with little influence on the pattern outside of this range. Through 
such adjustments it might be possible to obtain more refined estimates 
of atmospheric composition. However, no attempt to derive such refined 
estimates was made in the present study, in view of the inherent lack of 
precision involved in the determination of the transmission character- 
istics of the atmosphere in the present manner* 

It appears from these calculations, that the atmospheric inversion 
observed by Gillett, et al. in the 1306 cm”* CH 4 band should also be de- 
tectable from high-resolution measurements across the 3020 cm ^ band 

of CH, . 
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Microwave Spectrum of Jupiter . 

The brightness temperature of Jupiter at microwave wavelengths 
( 0,1 cm to 10.0 cm) was also computed for each thermal structure model. 
Strong atmospheric absorption occurs on Jupiter between 11,900 Me and 
40,000 Me in the microwave region due to the NH^ inversion effect. In 
these calculations , the ninety-two NH^ inversion lines which have been 
observed (Townes and Schawlow, 1955; Nishikawa and Shimoda, 1955) have 


been taken into account. 


At each frequency of interest, the molecular absorption aaeffi cr.t 
due to each line was computed following Townes and S chav: low (l nr > r ) t 
assuming the VanVleck - ’’eisskopf line share. The half-width was 
assumed to be due to foreign broadening of IT by K and lie . Accordin' 
to ^nderson (1950), the most important interaction in this case is be- 
tween the quadrupolar momentum of and the induced dipolar momentum 
of I?2 and He, and the cross-section for this interaction has been com- 
puted following his development. By sunn in a the contributions from each 
KHj line, as well as the contribution from nressure-induced absor- 
ption, the total opacity at each frequency was obtained# 

A/ith a knowledge of the total absorption coefficient at each point 
in the atmosphere, the brightness temperature was determined by standard 
radiative transfer techniques. The results for thermal Models 1, 2 and 
3 are shown in Figures 4 and 5, where a comparison with observations 
(tabulated by Law and Staelin, 1969, and Goodman, 1969) is made. Also 
included are the recent observations of Welch (1969). As can be seen 
from Figure 4, general agreement with the observations is obtained at 
wavelengths shorter than 3 cm in the 0.1 cm to 10.0 cm range. The 
minimum in brightness temperature observed between 1 cm and 2 cm in the 
region of strong NH^ absorption is well replicated in the three cases 
shown. At wavelengths longer than 3 cm, divergence between observed and 
computed brightness temperatures is found due to a non-thermal component 
in the observations which increases as the square of the wavelength. 

The computed thermal component increases with wavelength in a similar 
fashion due to the assymetric form of the assumed line shape. The 
VanVleck - ’.'eisskopf profile produces a slowly varying absorption coef- 
ficient at wavelengths distant from the line center on the short -wave- 
length wing, while on the long -wavelength wing the absorption coeffi- 


cient varies inversely as the square of the wavelength. 

As wavelength increases, the thermal component of brightness tem- 
perature computed here converges asymptotically to a value character- 
istic of the lower boundary temperature of the model, since sources of 
radiation below the top of the dense cloud layer were not considered. 
Therefore, in the present study, the thermal brightness temperatures 
near 10 cm wavelength have been underestimated. In fact, an observa- 
tional value of 260 •- 26K was deduced for the thermal component at 10. 4 
cm wavelength by Berge (1966). 

A more detailed depiction of brightness temperatures in the region 
of strong absorption is shown in Figure 5 for Models 1, 2 and 3 • Al- 
though the minimum temperatures computed at the center of the NH^ band 
agree well with the observations, the computed temperatures are warmer 
than those observed in the 0.95 - 1.05 cm and 1.35 - 1.45 cm ranges by 
some 10K. The reasons for this discrepancy are unclear. The most 
interesting aspect of this figure is the evidence of the atmospheric 
inversion which appears in the brightness temperature curves of Models 2 
and 3. In Model 1, where no inversion is present, the lowest brightness 
temperatures are found at wavelengths where strongest absorption occurs. 
In Models 2 and 3, however, the atmospheric inversion produces elevated 
temperatures (by as much as 5K) at these wavelengths. These models are 
consistent with the findings of Wrixon (1969) who reports evidence of a 
thermal inversion in recent microwave observations of Jupiter. 

Conclusions . 

These calculations indicate that thermal models of atmospheric 
structure, computed for boundary conditions suggested by recent obser- 
vations of Jupiter, give rise to infrared and microwave spectra which 


are in good agreement with measurements of t>e Jovian brightness 
temperatures in these spectral regions. Although the present method of 
evaluating atmospheric transmissions in t l e infrared rep-ion is probably 
not precise enough to permit a meaningful refinement of atmospheric 
composition estimates, it is felt that the models derived here are 
fundamentally sound, and sufficiently accurate to reveal several basic 
features of the Jovian atmosphere: 

There is an extensive "troposphere" overlying the dense clouds on 
Jupiter, in which the atmospheric structure is determined by dynamical 
processes, 

A thermal inversion of about 30K, produced by solar infrared heat- 
ing in the 3020 cm“^ CH^ band, is present in the mesosphere, as suggested 
by Gillett, et al, (1969) on the basis of their observations. 

The effective temperature of the atmosphere is significantly higher 
than that computed on the basis of solar eouilibrium. A model with a 
temperature and pressure at the cloud level of 220K and 4*6 atm, 
respectively, while producing the best overall agreement with t’*e 
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observed infrared brightness temperatures, is characterized by a radia- 
tive loss^ to space which is approximately twice that reciuired to 
balance the interacting solar flux for an albedo of 0 , 45 . An internal 
source of energy on Jupiter is therefore implied. 
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TABLE I 


MODEL CHARACTERISTICS 


(H 2 :He:CH 4 :NH 3 


5.0:1.0:0.005:0.001) 


odel 

H Pressure 
at Cloudtop 
(Atm) 

Cloudtop 

Temperature 

(°K) 

1 

2.0 

220 

2 

4.0 

230 

3 

4.0 

220 

4 

4.0 

210 

5 

6.0 

220 
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TABLE II 



AMOUNTS 

OF ATMOSPHERIC 

CONSTITUENTS 



ABOVE 

VARIOUS TEMFEEA 

TIIRS LEVELS 




MODEL 3 




Ho 

He 

CH, 

NH, 

T(°K) 

2 


4 

3 


(km-atm) 

( km-atm) 

(m-atm) 

(m-atm) 

220 

150 

30 

150 

21 

200 

114 

23 

114 

14 

ISO 

SO 

16 

SO 

7.6 

160 

55 

11 

55 

2.6 

140 

36 

7.2 

36 

.15 
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Figure Captions 


Figure 1. 
Figure 2, 

Figure 3. 
Figure 4* 


Figure 5. 


Thermal models for the Jovian atmosphere. 

Density distributions for H , He, CH. and NH_ corresponding 

2 4 3 

to thermal Model 3. 

Computed and observed infrared brightness temperatures of 
Jupiter, The dots correspond to the measurements of 
Gillett, et al, (1969), while the vertical bar at 20 
wavelength corresponds to t^e range observed by Low (19'$). 

Computed and observed microwave brightness temperatures of 
Jupiter, The computed values (indicated by solid curves) 
include only t’-e thermal component, while the observed 
brightness temperatures (indicated by vertical bars) include 
a non-thermal component which becomes ouite large at wave- 
lengths beyond 3 cm. 

Computed and observed microwave brightness temperatures of 
Jupiter near the center of the NH^ inversion band. Vertical 
lines at the top of the figure indicate the position of the 
strongest lines in this region. 
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